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' Abstract Wc analyze in situ measurements of solar wind velocity obtained by 

\ the Advanced Composition Explorer (ACE) spacecraft during the solar activity 

^ ' cycle 23. We calculated a robust complexity measure, the permutation entropy 

^1 ('S') of solar wind time series at different phases of a solar activity cycle. The 

' permutation entropy measure is first tested on the known dynamical data before 

its application to solar wind time series. It is observed that complexity of solar 
wind velocity fluctuations at 1 AU shows hysteresis phenomenon while following 
the ascending and descending phases of the activity cycle. This indicates the 
presence of multistability in the dynamics governing the solar wind velocity over 
a solar activity cycle. 

O ' 1. Introduction 

^ The corona changes its shape enormously during a solar activity cycle resulting 

in temporal and structural properties of the solar wind velocity variation in 
a solar cycle ( [Schwenn, 2007 ). In situ solar wind plasma observations show 
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^ ■ that its local properties at 1 AU are modulated by the solar activity cycle 

^sD ' ( [Hapgood et ai, 1991[[Rrchardson and Kasper, 2008[ ). Around solar activity min-[ 

imum, the structures of the corona and the solar wind are rather simple and re- 
I main so for several months ( [Schwenn, 2007 ). At solar activity maximum, slow so- 

. ■ lar wind dominates at all helio latitudes ( McComas, Gosling, and Skoug, 2000[ ) . 

(N ■ Recently it was reported that slow solar wind velocity just before the maximum 

of the solar activity cycle is least correlated to data obtained from the rest of 
the solar activity cycle ( |Suyal, Prasad, and Singh, 2011 1. 



Observations of solar wind velocity made by different spacecraft have been an- 
alyzed in considerable detail and reported by several authors. Values of complex- 
^ ■ ity measures such as entropy ( [Macek and Obojska[[1997.,1995l[Macek and Redaellil l 
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MTOllRedaelli and MacekllMJT]) . correlation dimension dMacek and Obojska, I997H 



Gupta et al, 2008 



and Lyapunov exponents ( |Macek and Obojska, 1998 Redaelli and Macek, 2001H 



Gupta et al., 2008) show that solar wind velocity fluctuations are a consequence 
of complex nonlinear dynamical processes. Inherent changes in the dynamics gov- 
erning the solar wind velocity at 0.3 AU have been observed ( [Gupta et al., 2008| . 
IMilano et al.\ (j2004p analyzed magnetic and bulk velocities measured by ACE 
and found an anisotropy in the velocity, magnetic, and cross helicity correlation 



functions and power spectra. [McComas, Gosling, and Skoug (|2000p analyzed 



Ulysses observations to demonstrate that the mid-latitude solar wind structure 
becomes increasingly complex as solar activity increases. [Consolini, Tozzi, and de Michelisj 
(|2009p investigated the emergence of spatio-temporal complexity in the 11-year 
solar cycle monitored by sunspot activity. They showed that spatio-temporal or 
dynamical complexity is an intrinsic property of the solar cycle. Using informa- 
tion entropy approach to the sunspot number time series, they showed how the 
dynamical complexity increases during the maximum phase of the solar cycle. 

Hysteresis occurs in several phenomena in physics, chemistry, biology, and en- 
gineering. It is a nonlinear phenomenon observed in systems from diverse areas of 
science, e.g., electromagnctism, electro-plasticity, superconductivity, and granu- 
lar motion (Bertotti, 1998,|Guyer, TenCate, and Johnson, 1999|[Katzgraber et al., 2002H 



Zharkov, 2002). This phenomenon occurs when a nonlinear system has at least 
two coexisting stable states in the hysteresis region where the system is found to 
depend on the history of the dynamics. Here, in dynamical systems, history cor- 
responds to the system's initial conditions. For example, hysteresis is observed in 
van der Pol system. Duffing system ( Thompson and Stewart, 1986D , and Lorenz 



system ([Alfsen and Fr(/)yland, 1985]) . The hysteresis has been observed in cou- 



pled nonlinear systems as well ( Prasad et al., 2005] ) 



Hysteresis phenomenon has been noticed in various solar indices. [Bachmann and White ] 
((T994)) observed the presence of hysteresis patterns among many pairs of activ- 
ity indices during solar cycle 21 and 22. They found that this hysteresis can 
be expressed approximately as a hierarchy of delay times behind the leading 
index, the sunspot number. [Jimenez-Reyes et aL|P998p analyzed the low-degree 
p-mode frequency shifts and solar activity indices (radio flux at 10.7 cm and 
magnetic index) over solar cycle 22 and observed a hysteresis phenomenon. 
IMoreno-Insertis and Solankil (|2000p suggested that high latitude fields are neces- 
sary to produce a significant difference in hysteresis between odd and even-degree 
p-modes frequencies. [Tripathy et al.\ (|2000p reported that the intermediate de- 
gree p-mode frequencies of solar cycle 22 show a hysteresis phenomenon with 
the magnetic in dices whereas no such effect exists for the radiative indices. 



Ozgiig and Atag (|200ip showed the presence of hysteresis between the solar 
flare index and some solar activity indicators such as total sunspot area, mean 
magnetic field, and coronal index during solar cycles 21 and 22. They found that 
these indices follow different paths for ascending and descending phases of the 
solar cycles while saturation effect exists at the extreme phases. 

In the present work, we attempt to understand the dynamics of ascending and 
descending phases of a solar cycle. We use permutation entropy {S) of hourly 
averaged solar wind velocity time series to capture the complexity trend over 
a cycle. We use the data obtained from ACE during 1998 — 2010. This period 
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belongs to solar activity cycle 23. We use permutation entropy (S) to detect the 
hysteresis in a dynamical system and calculate it for time series obtained from 
simulated as well as solar wind velocity data. 

In the next section, we review the algorithm to calculate the permutation 
entropy (S) of a time series. In Section |3] we describe how S detects the mul- 
tistability present in the modeled dynamical system. In Section |4] we analyze 
solar wind data using permutation entropy. This is followed by conclusions in 
Section [51 



2. Permutation Entropy 

Lyapunov exponent, entropy, and fractal dimension are well known and exten- 
sively used measures to detect dynamical changes in a time series obtained from 
a complex system ( |Kantz and Schreiber, 2004| ). Calculation of these quantities 
needs phase space reconstruction for which we need to know the parameters 
like embedding dimension and delay ( [Packard et al., 1980p . However, in prac- 
tice, it is very difficult to get accurate parameters, particularly for noisy data. 
Permutation entropy can be used to compare two or more time series and dis- 
tinguish regular, chaotic and random behavior for small and noisy time series 



(Bandt and Pompe, 2002). It quantifies not the only randomness but also the de- 



gree of correlational structures of a time series ([Rosso et ai, 2007). Permutation 



entropy is conceptually simple and computationally very fast ( Cao et al., 2004| ). 
It can be effectively used to detect dynamical changes in a complex system. The 
detailed algorithm to calculate permutation entropy is described below : 

Let us consider a time series Xi,i = 1,2, N and embed it in an m-dimension- 
al space ( [Packard et ai, 1980[ ). An embedded vector is written as 

For each embedded vector Xi, the m components can be arranged in an in- 
creasing order: [xi+ki < Xi+k2 < ■■■ < where kj can have any value 
from to m — 1. Hence, each vector Xi is uniquely mapped onto (fci, fc2, km) 
which is one of the m! permutations of m distinct symbols (0, l,...,m — 1). 
When each permutation is considered as a symbol, the reconstructed trajectory 
in the m-dimensional space is represented by a symbol sequence. We denote 
the probability distribution for the distinct symbol by Pi, P2, Pj , where j < 
ml. The permutation entropy for the time series Xi,i = 1,2,...,N is defined 
( [Bandt and Pompe, 2002[ ) as the Shannon entropy for the j distinct symbols 



Here S{m) attains the maximum value ln(77i!) when Pj = for uniformly 

distributed data. Therefore, normalized permutation entropy is written as 

in(m!j 
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Figure 1 . Schematic configuration of coupled oscillators represented by Equation lO . 

where the value of 5* lies between and 1. Smaller value of S indicates a 
more regular time series. If m is too small there are very few distinct states 
and this scheme will not work. A value of m = 5,6, or 7 is suitable to calcu- 
late the permutation entropy for detecting the dynamical changes in a system 
dBandt and Pompe, 20021 |Cao et al, 2004[ ). 

3. Dynamical Hysteresis 

Hysteresis occurs when a nonlinear system has at least two existing stable states 
in a hysteresis region. Out of multiple states, the system attains a state de- 
pending on the history of the dynamics, e.g., on the initial conditions of the 
system ( [Thompson and Stewart, 1986| [Prasad et ai, 200"5l ) . We consider one of 
the example from lPrasad et aL\ l|2005p , where the dynamical hysteresis in coupled 
oscillators have been observed in terms of stability of the system using Lyapunov 
exponents in a wide range of parameter space. It is conceptually difficult to 
estimate the Lyapunov exponent for a small and noisy data set. In order to 
see the hysteresis region using the permutation entropy, we consider Rossler- 
type ( [Gaspard and Nicolis, 1983| coupled oscillators, as shown schematically in 
Figure 1. The model equations for this system are 



dxj (t) 

dt 
dytjt) 

dt 
dzjit) 

dt 



= WiXi + ttiUi, 

= PiXi + Zi{xi - j/j), (3) 



where i,j=l,2,3,4 (number of oscillators^ 4), wi — 1.005, i/;2 = = W4 = 
0.995, ai = 0.38, /3i = 0.3, and ji = 4.5. At these set of parameters all the 
individual systems show chaotic oscillations ( [Prasad et ai, 20Q5[ ). The coupling 
functions are Fi = e{x2 + X3 + X4 — 3xi), F2 = e{xi — X2), F3 = e{xi —x — 3), and 
F4 ~ e(xi — X4) while e is the coupling parameter. In our numerical calculation, 
we used the fourth-order Runge-Kutta integrator with integration time step 0.01. 
[Prasad et al.\ (|2005p have shown the hysteresis behavior in the largest Lyapunov 
exponent with respect to the coupling parameter e, and predicted the coexistence 
of more than one state of different type of stability in the hysteresis region. In 
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Figure 3. The time series xi(t) in regions (a) Ai(e = 0.05), (b) A2(e = 0.20), (c) 
B2{e = 0.20), and (d) Bi{e = 0.05), respectively of Figurc[2] The right-hand-side panel shows 
the corresponding projection of trajectories into the x — y plane. 

order to check the presence of hysteresis and dynamical changes in terms of 
permutation entropy S, we consider the x-component of the time series of the 
first oscillator. We take embedding dimension m = 5 (although higher values of 
m gives similar results) and sampling time r — 300 (first minima of the auto- 
correlation function). 750 data points arc used to calculate S. Figure [2] shows a 
plot of S' as a function of increasing and decreasing coupling strengths e. The 
arrows show the corresponding increasing and decreasing coupling strengths. The 
hysteresis region is visible for different values of e. Figures [SKa — d) show the time 
series xi{t) and the projection of the first oscillator in the x ~ y plane at regions 
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corresponding to Ai{e = 0.05), ^2(6 = 0.2), ^2(6 = 0.2), and Bi(e = 0.05) of 
Figure [5] Although the trajectories look similar, the corresponding permutation 
entropy (S) shows that they are dynamically different. Importantly, the permu- 
tation entropy captures the hysteresis loop as shown with Lyapunov exponents 
in IPrasad et aD(|2005)) . 



4. Permutation Entropy of the Solar Wind Data 

We use hourly averaged solar wind velocity data at a distance of 1 AU ob- 
tained from the Solar Wind Ion Composition Spectrometer (SWICS) on ACE 



(http://www.srl.caltech.edu/ACE/ASC/level2/). This data set corresponds to 



the years 1998 to 2010. Since the data are not continuous we split the data 
set into 18 continuous time series to cover most of the solar cycle 23. The 
details of the time series used are given in Table 1. In addition. Table 1 (last 
column) contains the corresponding sunspot numbers for each time series taken 



from http://sidc.oma.be/sunspot-data/ Figure|l]shows the solar activity during 



activity cycle 23. The position of serial numbers corresponding to the solar wind 
data sets, is shown on the smoothed monthly averaged sunspot index curve. 

Table 1. Hourly averaged solar wind velocity data measured 
by the ACE spacecraft in the years 1998 to 2010: Initial time 
(Ti), Number of Data poin ts (A'^), and Sunspot numbers (SSN ). 
Sunspot numbers are from [http: / / sidc.oma.be/sunspot-data/] 



S.No. 






JV 


Sunspot number (SSN) 


1 


1998 


25 


1665 


54 


2 


1998 


71 


1314 


69 


3 


1998 


91 


789 


75 


4 


1999 


11 


1927 


84 


5 


1999 


41 


2716 


92 


6 


1999 


73 


1139 


104 


7 


1999 


94 


1404 


111 


8 


2000 


15 


3075 


117 


9 


2000 


88 


1580 


113 


10 


2001 


40 


1315 


110 


11 


2002 


03 


3329 


113 


12 


2002 


41 


3417 


107 


1.3 


2003 


03 


789 


81 


14 


2003 


24 


1840 


72 


1.5 


2003 


90 


2633 


54 


16 


2005 


33 


2541 


30 


17 


2006 


45 


2278 


16 


18 


2008 


39 


2108 


3.5 



We calculate the normalized permutation entropy S for all 18 time series 
using embedding dimension m ~ 5 and window length of 750 points. We use six 
overlapping windows and then calculated the average value of S. The standard 
deviation in S is given by 



n 
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Figure 4. Smoothed monthly sunspot number from 1996 to 2010. Serial numbers 1 — 18 
correspond to the solar wind data sets given in Table 1. 




where n is the number of overlapping windows, S'avg is the average permutation 
entropy of n windows, and a estimates the statistical uncertainties in the average. 
Figure [5] shows the plot of S'avg of solar wind time series vs. sunspot number 
corresponding to the time of solar wind observation. The error bars are estimates 
of the statistical uncertainties in these averages. 

It is noticeable that the path followed by the permutation entropy of the solar 
wind data in the increasing phase of the solar activity cycle is different from that 
of the decreasing phase, i.e., a hysteresis phenomenon is present. It shows that 
although sunspot data set is almost symmetric around the peak of the solar 
activity cycle yet permutation entropy of solar wind data follows different paths 
in increasing and decreasing phases of the solar activity cycle. 
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5. Conclusions 

The long term sunspot time series shows an average cycle length of 11 years. 
Although nearly periodic, the period as well as amplitude of the cycle varies 
irregularly. Apart from the sunspot data, the intrinsic irregularity of the solar 
cycle is seen in other observable variables like surface flows, solar irradiance (solar 
constant), the solar wind, and so on. The Maunder Minimum and several an- 
cient periods from solar proxy-data suggest that the Sun exhibits quasi-periodic 
or intermittent behavior ( |Feminella and Storini, 1997[ ). Owing to changes in 
magnetic activity, many aspects of the solar wind change over a solar cycle, 
including the speed, the density, the dynamic pressure, the composition, and the 
temperature ( [Richardson and Kasper, 2008[ ). 

In this paper, we present the analysis of solar wind velocity data during 
the solar activity cycle 23. We obtained 18 different time series of hourly av- 
eraged solar wind velocity, measured by ACE spacecraft at 1 A.U. To quantify 
the randomness of these time series, we use a robust, conceptually simple and 
computationally efficient measure called permutation entropy. A smaller value of 
permutation entropy indicates a more regular time series. We observe that as the 
solar cycle 23 progresses towards maximum, the permutation entropy increases, 
saturates around the peak of activity and then decreases as the activity of cycle 
23 subsides. We also note (Figure. [5]) that the value of permutation entropy 
follows different paths in the ascending and descending phases of the solar ac- 
tivity cycle. In addition, while the ascent is fluctuating, the descent is smooth. 
This hysteresis phenomenon shows the multistability in the dynamics of solar 
wind, over the solar activity cycle. The behavior is similar to the one observed 
for hysteresis phenomenon of other solar indices ([Bachmann and White, 1994 



Jimenez-Reyes et al, 1998 Ozgiig and Atag, 2001) and confirms (c/. Consolini, Tozzi, and de Michelis[ [ 



2009P that spatio-temporal or dynamical complexity is an intrinsic property of 



the solar cycle. 

Acknowledgements The authors thank the ACE Science Center and instrument teams for 
making available the ACE data used here. VS and AP thank CSIR for SRF and DST Govt, 
of India for financial supports respectively. 



References 

Alfsen, K.H., Frc/iyland, J.: 1985, Phys. Scr. 31, 15. 
Bachmann, K.T., White, O.R.: 1994, Solar Phys. 150, 347. 
Bandt, C, Pompe, B.: 2002, Phys. Rev. Lett. 88, 174101. 

Bertotti, G.: 1998, Hysteresis in Magnetism: For Physicists, Materials Scientists, and 

Engineers, Academic Press, New York. 31. 
Cao, Y., Tung, W., Gao, J.B., Protopopcscu, V.A., Hively, L.M.: 2004, Phys. Rev. E 70, 

046217. 

Consolini, G., Tozzi, R., de Michelis, P.: 2009, Astron. Astrophys. 506, 1381. 
FemincUa, F., Storini, M.: 1997, Astron. Astrophys. 322, 311. 
Gaspard, P., Nicolis, G.: 1983, J. Stat. Phys. 31, 499. 

Gupta, K., Prasad, A., Saikia, I., Singh, H.P.: 2008, Planet. Space Sci. 56, 550. 
Guyer, R.A., TenCate, J., Johnson, P.: 1999, Phys. Rev. Lett. 82, 3280. 

Hapgood, M.A., Lockwood, M., Bowe, G.A., Wills, D.M., Tulunay, Y.K.: 1991, Planet. Space 
Sci. 39, 411. 



SOLA: vinita_et_al.tex; 23 December 2011; 2:19; p. 8 



Hysteresis in a Solar Activity Cycle 



Jimenez-Reyes, S.J., Rcgulo, C, Palle, P.L., Roca Cortes, T.: 1998, Astron. Astrophys, 329, 
1119. 

Kantz, H., Schreiber, T.: 2004, Nonlinear Time Series Analysis, Cambridge University Press, 
Cambridge, 10. 

Katzgraber, H.G., Pazmandi, F., Pike, C.R., Liu, K., Scalettar, R.T., Verosub, K.L., Zimanyi, 

G.T.: 2002, Phys. Rev. Lett. 89, 257202. 
Macek, W.M., Obojska, L.: 1997, Chaos Solitons Fractals 8, 1601. 
Macek, W.M.: 1998, Physica D: Nonlinear Phenomena 122, 254. 
Macek, W.M., Obojska, L.: 1998, Chaos Solitons Fractals 9, 221. 
Macek, W.M., Redaelli, S.: 2000, Phys. Rev. E 62, 6496. 

McComas, D.J., Gosling, J.T., Skoug, R.M.: 2000, Geophys. Res. Lett. 27, 2437. 

Milano, L.J., Dasso, S., Matthaeus, W.H., Smith, C.W.: 2004, Phys. Rev. Lett. 93, 155005. 

Moreno-Insertis, F., Solanki, S.K.: 2000, Mon. Not. Roy. Astron. Soc. 313, 411. 

Ozgiig, A., Atag, T.: 2001, In: Brekke, P., Fleck, B., Gurman, J. B. (eds.). Recent Insights into 

the Physics of the Sun and Heliosphere: Highlights from SOHO and Other Space Missions, 

lAU Symp. 203, 125. 

Packard, N.H., Crutchfield, J. P., Farmer, J.D., Shaw, R.S.: 1980, Phys. Rev. Lett. 45, 712. 
Prasad, A., lasemidis, L.D., Sabesan, S., Tsakalis, K.: 2005, Pramana 64, 513. 
RedaeUi, S., Macek, W.M.: 2001, Planet. Space Sci. 49, 1211. 
Richardson, J.D., Kasper, J.C.: 2008, J. Atmos. Solar-Terr. Phys. 70, 219. 
Rosso, O.A., Larrondo, H.A., Martin, M.T., Plastino, A., Fuentes, M.A.: 2007, Phys. Rev. 
Lett. 99, 154102. 

Schwenn, R.: 2007, In: Baker, D. N., Klecker, B., Schwartz, S. J., Schwcnn, R., von Steiger, R. 
(eds.). Solar Dynamics and its Effects on the Heliosphere and Earth, Springer, New York, 
51. 

Suyal, v., Prasad, A., Singh, H.P.: 2011, Submitted to Planet Space Sci.. 

Thompson, J.M.T., Stewart, H.B.: 1986, Nonlinear Dynamics and Chaos, John Wiley and 
Sons, New York. 123. 

Tripathy, S., Kumar, B., Jain, K., Bhatnagar, A.: 2000, J. Astrophys. Astron. 21, 357. 
Zharkov, G.: 2002, J. Exp. Theor. Phys. 95, 517. 



SOLA: vinita_et_al.tex; 23 December 2011; 2:19; p. 9 



SOLA: vinita_et_al.tex; 23 December 2011; 2:19; p. 10 



